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Abstract 



We study correlation functions of parafermionic currents and disorder 
fields in the Zn symmetric conformal field theory perturbed by the first 
thermal operator. Following the ideas of Al. Zamolodchikov, we develop 
'■ for the correlation functions the conformal perturbation theory at small 

\ scales and the form factors spectral decomposition at large ones. For all 

N there is an agreement between the data at the intermediate distances. 
We consider the problems arising in the description of the space of scaling 
fields in perturbed models, such as null vector relations, equations of mo- 
tion and a consistent treatment of fields related by a resonance condition. 



> 

■ 1 Introduction 

, The calculation of correlation functions is one of the most interesting problems in 

' the two dimensional integrable quantum field theory. Complete exact solutions 

, of this problem were found for free field models, as, for example, Ising model, 

and for the Conformal Field Theories (CFT). 

In massive integrable field theories an elegant way of studying behaviors of 
two point correlation functions was proposed by Al. Zamolodchikov in Ref. pQ. 
In the case of the perturbed Lee- Yang model he studied short and long dis- 
tance asymptotics of the two point correlation function of primary fields. The 
short distance behaviors of correlation function were investigated by developing 
the infra red safe perturbation theory [1] , based on the knowledge of the exact 
Vacuum Expectation Values (VEVs) of local fields found by non-perturbative 
methods [21 [31 H] . The correlators in the infra red region were given by using a 
form factor spectral decomposition [5l [6] . A very good agreement between the 
asymptotics at the intermediate distances Mr ~ 1 allowed to claim that cor- 
relation functions in Lee- Yang Model are effectively described in this approach 
at all distances. A further development has shown that the proposed in Ref. 
[T] method is a simple and effective tool for an analysis of basic properties of 
correlation functions in different integrable massive models. 
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In the present paper we continue our studying [7J IB] of scaling fields in the 
parafermionic CFT [9j [10] with the central charge 



c=2(N-l)/(N + 2), N = 2,3,---, 



(1) 



perturbed by the first thermal operator Si 



A = Ac ft + A / d 2 x £\(x) 



(2) 



The resulting theory is integrable and Z/v symmetric. Depending on the sign of 
A the system is in the ordered or disordered phase. We fix A > phase where the 
Zm symmetry is destroyed and vacuum expectation values of disorder operators 
are non-zero. In this scaling model 112) . we study, within the approach 
of Ref. PP, the correlation functions of disorder operators and parafermionic 
currents. These objects, as well as correlators of some VF-algebra descendants, 
are not so easy for a direct investigation, mainly because of the resonances, 
which appear in the construction. Namely, following the procedure of Ref. pQ , 
we need to study the situations, where the scaling dimensions D a and Db of 
some fields O a and Ob satisfy the condition 



In this case, the field O a has the n-th order resonance with the field Ob and 
there is an ambiguity in defining the renormalized field O a 



This typically results in a logarithmic scaling of the field O a [!]■ We observed 
that according to our general formulae the form factors of fields satisfying ([3]) 
formally coincide as functions of rapidities. This creates a problem of defining 
form factors of such scaling fields, since it is expected from the general settings 
[T3J[T2], that there is a one-to-one correspondence between conformal and scaling 
fields. We suppose to discuss a general prescription for form factors of the fields 
possessing resonances in a separate publication. In this paper we consider some 
examples of the phenomena. 

From the other side, studying short distance behaviors of the correlators, we 
found, that the terms in the perturbative expansion, which correspond to fields 
with the condition ([3]), formally diverge for integer parameters N. To proceed 
further with these cases, we use the fact, that our exact expressions are defined 
for the models with arbitrary number N . We provide an analytic continuation 
over this parameter and obtain finite results for the correlators. We check, that 
our prescription leads to correct expressions for the Ising model (N — 2 case) 
[T3] , as well as for other known cases. 

A matching between the long and short distance asymptotics, which we 
found for all N, can also be considered as an additional confirmation, support- 
ing the consistency of the proposed form factors, the expressions for the short 
distance expansions, for VEVs and for the normalizations of the scaling fields. 



D a = D b + 2n(l- A £l ), n>0. 



O a ^O a + const A™ O b • 



(4) 
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We choose for studying the Zjy invariant Ising model, since it has several 
nice properties. It is related with a statistical system, which have a simple 
and clear description [5J [TB] and many important physical applications, while 
the operators of this lattice model, in general, have an interesting quasi-locality 
property [T7J [TH] in the scaling limit. We recall the notions of the lattice Zjv 
model in the section 2. In the critical points the lattice model is described by the 
parafermionic CFT with the central charge ([T]). In section 3, we introduce neces- 
sary definitions and collect basic facts about this CFT, concentrating attention 
on the algebraic structure of the space of conformal fields. This information is 
essential for further analysis of the space of form factors of local fields in the case 
of massive integrable model, which we consider in the section 4 as the model of 
QFT, defined as scattering theory with the simple Zm symmetric S matrix [T2] . 
We discuss an algebraic prescription for the form factors of the scaling fields. 
In particular, we consider the actions of deformed parafermionic currents in the 
space of form factors and study related questions, such as deformation of the 
quantum equations of motion, null vector relations and a prescription for form 
factors of fields, satisfying the resonance condition ([3]). Finally, the correla- 
tion functions of scaling fields of the massive integrable model ([2]) are discussed 
in the section 5. We concentrate our attention to the conformal perturbation 
theory and demonstrate, how to apply the quantum equation of motion to the 
computing coefficients of the perturbation theory. We discuss the regularization 
prescription for the resonances, appearing in the perturbation theory and also 
provide results of numerical computations. 



2 Lattice Z^-Ising model 



In this section we recall following Ref . [5] basic definitions of the two dimensional 
lattice model, generalizing the well-known ^-symmetric Ising model to the case 
with the Zjy, {N = 2,3,...} symmetry. 

We consider the model of the statistical mechanics defined at the square 
lattice. Let spin variables a are associated with the cites of the lattice and take 
values in the group Zn, i.e. a = 



(k = 0,...,N-l), where u = exp (^p) 



Fig 1. Zm Ising model at the lattice. 

We study the local theory where Boltzmann weig ht e- n ^-' J '"> = W(a, a') de- 
pends on spins a and a' situated at neighboring sites. The partition sum of the 
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model is, by definition 

z = e n viw). (5) 

spins edges 

In Z/v Ising model Boltzmann weights satisfy the Zn symmetry W(u, <r') = 
W(o;<7, wcr'), and the reality condition W(a,a' means, that, 

up to a normalization constant, the function W(<7, a') has the form 

JV-l 

W(a,a') = m(^v')\ Wo = 1 , (6) 

where the parameters W& of the model are real non-negative numbers satisfying 
the equation Wfc = Wjv-fe- 

An important information about the model is encoded in a set of its correla- 
tion functions. For example, the correlation functions of the spin operators are 
defined at the lattice as 

(<J kl (xi) ■ ■ ■ a ks (x s )) = i ^(xi)--- a k ° (x s ) [] W(a, a') . 

spins edges 

We will study a system in the thermodynamic limit, assuming appropriate peri- 
odic conditions at the infinity. The one point correlation functions (<7fc) serve as 
a measure of order at the system. Another important lattice operators are the 
disorder operators [TUl [TT1 [18] . Consider a directed path T going through points 
of a dual lattice and intersecting the bonds of the original lattice. Let all weights 
W s at the bonds, crossed by path T, are changed to become wi fc ' ) = W s uj sk . 
The presence of the dislocation along the path introduces in the system a frac- 
tional domain wall, that favors a discontinuity in the value of the neighboring 
spins by k. The partition function Zj? on the inhomogeneous lattice now differs 
from the original one. We interpret the dislocation as an insertion of two oper- 
ators /ifc and fjl = fiN-k situated at the sites of the dual lattice as it is shown 
in Fig.2. 
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Fig 2. Two point correlator of disorder operators fik and fil. 
By definition, the two point correlation function of conjugate disorder operators 
is given as 

(MfcOi)/4(^2)) r = ^ ^ r ' 
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Here 2:1,2 are the coordinates of the corresponding operators. This interpreta- 
tion turns out to be useful since the dependence on the contour T is not very 
strong. As in the Ising model [17] . one can make contour deformations by closed 
paths, without changing the correlation function. Applying this freedom, we use 
prescription of attaching contours from infinity to each of the points Xj of the 
dual lattice, as it is shown in Fig. 3. Now, the definition of two point correla- 
tion functions can be immediately extended to the multi-point case, including 
disorder operators with a total non-zero Zn charge. 
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Fig 3. Correlation function (/ifc /Ujv-fc}- 

In a general case of correlators of spin and disorder operators, only an absolute 
value of the correlation function remains to be a path independent, and one 
should fix relative position of contours [10] . This is because, the correlation 
functions of the product of spin and disorder operators, before and after a com- 
plete counterclockwise rotation of the disorder variable around order variable, 
differ by a phase uj~ m . Equivalently, this happens, whenever a ^-path crosses 
a (Ji variable. Operators with that properties are called mutually quasi-local 
with the exponent 7^ = —kl/N. 

Order and disorder parameters are basic operators in the theory. Other 
operators are constructed as their operator products. In general, if the distance 
\x — y\ between operators is much smaller than the correlation length, then the 
local operators ipi,k-, appearing at the operator product 

ai(x)Hk(y) = Ci k {x - y)ipi,k(x) H , 

do obey the parastatistics [18]. In particular, further we will study the corre- 
lation functions of two parafermionic currents if> = "01, 1 and ip^ in the scaling 
limit. 

Under the Krammers Wannier duality [17] order parameters become disorder 
parameters and wiseverse. The hyperplane of self-duality has the dimension 
[■j] . For the Ising model (TV = 2 case) and three state Potts model (TV = 3 
case), the system has a second order phase transition at the self-dual point. For 
TV = 4 it coincides with the well-known Ashkin- Teller model. For prime TV > 5, 
the Z N theory has, as a rule, three phases (ordered, disordered and Kosterlitz- 
Thoulcss phases). For non-prime TV the phase structure is more complicated. 
As an example, we consider the phase diagram of the Z§ model [16] . 
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Fig 4. The phase diagram of the Z5 model 

Here the phase I, the ordered phase, is characterized by the conditions (a{) 7^ 
and (fik) = 0. In the disordered phase II the situation is inverse (ct/) ^ and 
(lik) = 0. Finally, in the Kosterlitz-Thouless phase III expectation values of 
operators of both type are zero (cr;) = (^k) = 0. The line FB denotes the self- 
duality region. It contains two symmetrically situated "bifurcation points" C 
and C*. Along the line CC* the model have a first order phase transition and 
ordered and disordered phases may coexist. The points C and C* are integrable 
and critical |28[ 116] . The theory in these points has Zn x Zn symmetry and 
is described by the continuous parafermionic CFT constructed in Ref. [9]. In 
the scaling limit in the vicinity of a critical point the order 07 and disorder fj.^ 
operators, as well as parafermions ip, are described by the fields, depending on 
the continuous space parameters. We preserve for that fields the same notations 
as on the lattice. 

Due to the conformal invariance and infinite-dimensional symmetry of the 
critical theory [15], the analysis of its correlation functions, as well as the struc- 
ture of its space of states, simplify drastically. This was described in details 
in Refs. [9_1 [10]. Basing on CFT results, one can study basic properties of 
the correlation functions of the Zn models in a vicinity of the critical point 
by application of the conformal perturbation theory. Again, the easiest for an 
analysis perturbations are those, which are integrable, due to a presence of an 
infinite set of integrals of motion [13] . Different integrable perturbations of the 
conformal field theories [9J were studied in Ref. |11) . 

In this paper, we study a vicinity of the critical point C in the phase II, where 
the temperature deformation leads to an appearance of the finite correlation 
length and to non-zero vacuum expectation value of the disorder parameters [111 
[121 [20]. In the QFT language, that is the massive perturbation of parafermionic 
CFT ([2]) by the most relevant first thermal operator e\, which destroys the dual 
Zn symmetry and preserves the symmetry Zn . 

One of the questions, which we would like to address in our study, is a 
structure of the space of scaling fields. We expect, that the space of composite 
quasi-local fields constructed from operators <Jk and \i\ in the vicinity of critical 
point, will be essentially the same, as in the CFT. Let us note, that this state- 
ment was supported by several results on counting of local operators in various 
integrable models [HI [22] . To understand this problem, we tried to apply in 
Ref. [8] a knowledge on the algebraic structure of lattice operators, which is 
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based on a deformation of conformal algebras [23J 123 . According to this 

idea we recall first the structure of quasi-local fields in the CFT point. Then, 
we will try to apply the clear and simple algebraic scheme for the investigations 
of matrix elements of scaling fields in the basis of asymptotic states and also for 
studying their correlation functions off-criticality. 



3 Space of states in parafermionic CFT 

In this section, we recall basic facts [9J [TU] about the conformal field theory with 
the parafermionic symmetry which describes critical points of Zn Ising model 
[5] and related models [52] . 



3.1 Parafermionic symmetry 

In the conformal limit the order o k and disorder fj, k parameters, which determine 
long range correlations of spins and dual spins, have the anomalous dimensions 

Mk - N(N + 2)- (?) 

Their Zn and dual Zn charges, respectively, are equal to k. Under the action 
of the Zn symmetry spin fields transform as 

a k ^u kn a k , neZ. (8) 

The transformation law for operators fik under the action of the dual Zn sym- 
metry looks similarly 

n'eZ. (9) 

The spin and disorder fields are the basic operators in the theory. All other 
fields are constructed from them. The composite fields are naturally separated 
into families with the fixed additive Zn x Zn charges (k, I). The members of 
the family behave under Zn x Zn transformation as 

$ _^ W fcn-Mn' $j n ,n'eZ. 

The parafermionic fields ip k an d ip k generalizing the usual Ising model fermions 
appear in the OPE of the order and disorder fields 

<r k (z, zK(0,0) = \z\- id *z A *[ MO) + ■••], 
a k (z,z)nl(0,0) = \z\^z A nM0) + ■■■]• (10) 

These currents are holomorphic and generate the infinitely dimensional symme- 
try due to the conservation laws 

d^fc = 0, d z $ k = 0. (II) 
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We concentrate our attention to the simplest solution of the associativity con- 
dition for the operator algebra of currents, which corresponds to CFT with the 
central charge ([T]), and the conformal dimensions of currents 

M) 

The fields ip = tpi (and respective antichiral currents) are the basic ones in the 
parafermionic algebra. It is convenient for us to consider, as well, conjugate 
currents ipN-i = ip- 

In the conformal model the space of states splits naturally into a direct sum 
of subspaces with the specified Zn x Zm charge (k, I) 

{F} = ($ {F} [m - } , N>m,m>l-N, (13) 

where [m,m] = [k + I, k — I] , m + rft G 2Z. In these notations parafermionic 
currents and order-disorder fields belong to the following subspaces 

4> e {^} [2 ,0] - ^ g i F }[-2,0] . ( 14 ) 
^G{F} t0j2] , ^ € {F} [0 ,- 2 ] . 

°> e {F} [k k] , n k e {F} [fe _ k] . 

Conformal fields are classified according to the representations of the parafer- 
mionic algebra. The action of the parafermionic generators A v (A\) is defined 
by the OPE 

•4>{z)$[m,m] = Yl z ~^ +n ~ lA l±p--n $ [m,m] , (15) 

Notice, that, if &[ m ,m] € -F] m ,m] has the conformal dimensions (d, d) then the 
conformal dimensions of fields 

[m,m] 

G{F} 

[m+2,m] ; [m,m] 

G{F} 

[m— 2,m] ; 

are (d — v,d). 

Order and disorder fields are the primaries of the parafermionic algebra. For 
instance, the following equations hold for n > 

A^ +n Hk = A\_, +n+l Hk = 0, (16) 

All other fields of the model are obtained by the action of the currents ip, ip on 
the fields The space of states of the CFT decomposes into a direct sum of 
irreducible representations of the parafermionic algebra 

{F} = e^o 1 \»k] A ,A ■ (17) 

This space can be also obtained by an application of parafermionic generators 
to the order parameters Cfc. Moreover, in some situations we will need another 
infinite symmetry description of parafermionic CFT, namely, the W symmetry, 
which we recall in the next subsection. 
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3.2 W algebra symmetry 

The space of fields in parafermionic CFT allows a classification with respect to 
another infinite dimensional symmetry algebra, the so-called Wjy algebra [27] . 
The generators of the latter W 2 (z), W 3 (z), . . . appear at the operator product 
of parafermionic currents 



'(^3(0) + ^±HflW 2 (0)) + •••). (18) 



Here W2 currents with spin 2 generate Virasoro algebra with the central charge 
given by eq. (JTJ. The currents W3 have spin 3. They generate the whole algebra, 
including the higher spins currents, which are omitted in eq. (fT5]l . 

Currents of W algebra, and respective anti chiral currents, have zeroth Zjy 
charges. Acting on the highest weight fields, they create an irreducible repre- 
sentation. From the viewpoint of the Wn x Wn symmetry, each of the spaces 
\fik\A A expands into a direct sum of representations. Namely, let us denote as 
0- ' )Vfe the field with the minimal conformal dimension, that can be obtained 
by the I— times application of the parafermionic generators ip^ to /i£ 

(y^'/ifc = Ai- 3 - k ■ ■ ■ A- k Mfc ■ (19) 

N N N 

Its conformal dimensions are easily computed from (JTJ. Then, up to a normal- 
ization, the following relations take place 

*fc-ai,-fc+ar = ^ t ) , (^)W. U = 0,l,...,fc, (20) 
®l%i,-k-2i = W(? f ) W , 1,1 = 0, 1, N-k. 

(k) 

These fields &m'm are the W algebra primaries. The action of generators of 
W algebra on it create the irreducible representation [3?^ mWw- ^he explicit 
values of the conformal dimensions (dm , dm ) of these fields are given as follows 

,(k) (fc + l) 2 — 1 m 2 ^ , ^ . , 

d = ^— - — i , -m<k<m. (21) 

For the cases, when \m\ > k, we use the relation d^m = d^_^\ which follows 
from the Zjq symmetry condition 

- m-N,-N+fh ■ 



*(*) _ $ (oo) 



(2k) 

For example, the physically important energy fields ek = 3>o are among these 
primaries. Operators Sk are local with respect to all fields and have the confor- 
mal dimensions 

D k =k(k+l)/(N + 2). (23) 
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We refer to the paper [§] for further details. Let us only comment here, that the 
Z/v symmetric parafcrmionic CFT can be equivalently considered as a particu- 
lar case of the conformal field theory WAi p N , introduced in Ref. [37]. Namely, 
for the fixed value of N it has the smallest central charge (UJ among all ratio- 
nal unitary minimal conformal theories with the extended Wn algebra, which 
corresponds to the parameter p = N + 1 . 

We will use the conformal data, represented above in the next sections, to 
describe the perturbed conformal operators, for which we again will preserve 
the same CFT notations. 



4 Form factor approach 

Now we turn to the operators in the corresponding massive integrable theory @, 
which allows several equivalent descriptions. In the given section, we describe it 
as a two dimensional QFT model with the factorized scattering of Zm charged 
particles. The particles a <G {1, . . . , N~ 1} in the Zm (N — 2, 3, 4, . . .) symmetric 
models have masses [20, 

. r sin(7ra/iV) . 
M a = M — V / . (24) 
sin(7r/A0 v > 

The antiparticle is, by definition, identified with the particle N — a. The 
scattering matrix of the lightest particles a = 1 has a simple form |12) 



sinh(f + f) 
sinh(| - f ) 



*n(/3)= ■ (25) 



The S matrices for higher particles are also diagonal and can be extracted from 
Sn, according to a standard bootstrap prescription. For example, the scattering 
matrix between the particle 1 and the antiparticle is S' 11 t(/3) = Sii(i7r — /3). 

The knowledge of the exact spectrum (|24|) and the scattering matrix (|25|) 
allows to study correlation functions of the theory, by using its spectral decom- 
position into the series of form factors. The form factors 

(0(x)\p 1 ,...,p n ) ai _ tan , (26) 

of the scaling field 0(x) are matrix elements of this operator in a basis of 
asymptotic states, formed by particle creation operators. We assume, that the 
particles, labeled by a%, ■ ■ ■ , a„, have the rapidities • • • , f3 n . The functions 
(|2"6"|) should satisfy to some analytical properties and, also, to a set of functional 
equations, the so-called form factor axioms [3|6], to guarantee the (quasi)locality 
of scaling fields. 

An usual problem in the form factor approach (see, for example (29l [30l [31] 
for Zn Ising models case) is that it is not so easy to determine, which scaling field 
is described by the given functions, satisfying proper functional equations and 
analyticity conditions. Moreover, for a matching with short distance formulae, 
it is necessary to determine the normalization of the scaling fields, described by 
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form factors. In our construction |S], we proposed some algebraic approach 
to solve these problems. Still, there are many subtle questions in this direction. 
We would like to discuss some of them in the next sections. 

4.1 Free fields construction for form factors 

In Refs. [3 [8] we have followed the algebraic approach [33l [25J [23l [26l [34] to the 
form factors of scaling limit of the ABF model 32 . This lattice model falls into 
the same universality class, as the Zn Ising model. In the corner transfer matrix 
approach its hidden symmetry is a deformation of CFT symmetry algebras (J7J)- 
(j23|) . We explore algebraic maps in the space of form factors to produce exact 
expressions for form factors of scaling fields. Our basic prescription, derived 
in Ref. [8] for form factors, can be shortly re-formulated as following. We 
introduce the notations 

V ZSm N 0= ± 

B{P')t] m = -"^= £ b ewH-i-^^(^) (27) 

The explicit expressions for the form factors of fields (ijj') l {ip) l Hk in the per- 
turbed theory are given for m — k — 21, m — k — 21 as 

n n 

{{^) i m^\{P}AP'}) M = C&«IlB()8i)&nB+(^)^)> .(28) 

i i 

We assume in this equation, that the form factor is zero, unless the Zjy neutrality 
condition 2(n' — n) = m + fh is satisfied. The label n stands for the number of 
lightest antiparticles U, carrying the Zn charge -2, and n' means the number 
of particles 1 with the Z^ charge 2. The constant C,^^ is determined by 
the normalization of the scaling field. In what follows we will specify its value 
for spin-less fields and for parafermionic currents to be in agreement with the 
conformal normalization. In eq. (|28p . we used a shorthand notation for the 
state with these numbers of particles and antiparticles 



{/?'})(„,„,, = |/3i,...,/3 n X---X)it...it,i...i. (29) 

The symbol of ordered product of particle creation operators in eq. (|28|) is used 
for the object 

n 

Y[B(f3 ] )=B(p 1 )---B(l3 n ), (30) 

1 Matrix elements of operators Z a , Z\ are defined by the rules (31)-(32) below. A more 
explicit definition of these operators can be found in Ref. 1231 . 



11 



which is a linear combination of the products of exponential free bosonic fields 
Z±(f3), see R- ei - [22] for details. The expectation value of a product of operators 
Z a {P) and Z\{f3) over the Fock vacuum can be computed by applying the Wick 
theorem 

({z ai (A) • • • z an (p n )zl (/?;)••■ z\ m (/?;„))) (31) 

= n iz ai {Pi)z aj (&)» n $k w x 

n«^(fto<(A'))>, 
•J 

The contraction rules of two operators in this equation are determined in terms 
of the meromorphic functions C(/3), given in the appendix B, as following 

(we assume, that j3 — (5\ — fa) 

({zmzM)) = ((zUMzUM)) (32) 

sinh ^ 

= C t (/3)cosh(^-^(a + 6)). 
4.2 Parafermionic currents action 

Let us shortly comment the algebraic structures, encoded in the form factors 
prescription described above. We obtained the equation ([28]) starting from the 
thermal operators operator form factors case m — —m and n' = n |23j as a 
result of the parafermionic currents action [7] on the bosonic operators 

ZliJS) ^ Zl(p)e&-%r a , Z b ((3) Z 6 (/3)e"^- ^ b , 

Zt(P) n Zt(P)e-^° , Z b (P) n Z b (P)e^ b , 

ZUfi) ^ Zl((3)e-& + ^ a , Z b {0) ^ Z b (/3)e% + & b , 

&JJS) Zt(l3)e&-™ a , Z b ((3) ^ Z b ((3)e-^-^ b . (33) 

Such prescription determines a set of maps in the space of multi-particle form 
factors, for example H 

<O|*&|{0M/8 , }>(n,n') ^ (0|W<Ul{/3},{/3'}>(„-i,„') • ( 34 ) 

This relation, and similar relations for other parafermionic current actions, are 
understood as following. We assume, that the mult i- particle form factors for 

2 We use in eq. d 341) the conformal notations (VO^m'mi s i ncc this equation is a perturbed 
analogue of eq. (| 1 9 1> . 
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(k) 

the field ^ m \m are given as matrix elements of a linear combination of products 
of n operators Zb(f3) and n' operators Z\(f3). Then, the form factors of the 
field (V>)$mm be 

given by the modified bosonization prescription, where 
the number of antiparticles is one less, and the change (|33|) is provided for 
each of Z a and Z\ operators. One can see, that this prescription obviously 
agree with the equation ([28]) . The eq. (|33|) was obtained from the deformed 
parafermionic action. For this reason, the following identification was proposed 
for Z,[=0,l,...,fcmRef. 

(°l $ k-«,-*+2rl^>'^'»(».n') = (0\(^Y(W^m,{P'})(n,n') ■ (35) 

The validity of this equation for multi-particle form factors with m = — m was 
checked by comparing a prescription for the form factors (|27[) - (133l) with the form 
factors of the scaling fields in the left hand side of this equation, obtained in 
the prescription for the deformed W algebra primaries [33] . We also get a clear 
evidence on the correctness of this identification by studying the correlation 
functions of the order and disorder fields in Ref. [7 |f| We consider the equation 
(j35|) as an off critical analogue of the equation ([20)) . 

It is possible to continue further a studying of the structure of the space 
of form factors of scaling fields in an analogy with the correspondent algebraic 
description of CFT. We note here, that, together with the equation (|20|) . the 
form factor prescription ([27]) -([33 ]) also satisfy the charge conjugation condition 



(|2"2"|) . Using the Zn symmetry condition ip — < f > ^2 2 -JV' ^ * s P oss ible to derive 
the form factors of the parafermionic currents. A simple observation is that 

(2) 

these matrix elements are related with the form factors of the field $ 2 q = (4>) £ i 
as following (see also |29j ) 

yQl' ° -X' ■'.) (*aiw.{^» (n _ llB) . 

(36) 

This equation is an off critical analogue of the quantum equation of motion 
for the perturbed parafermionic currents 

J^M) = A^MeiOM) = \\f^<5>%{zrz) ■ (37) 

In the general multi particle case this equation follows from the trigonometric 
functions identities. Further we will use eq. (|3"T)) in the short distance expansion. 

Our next comment on a mapping between form factors and interpreting our 
form factors as matrix elements of the scaling fields is as following. Formally 
applying the prescription of Ref. |7 , we arrive to fact, that the null vectors 
equation (TTB]) survives after the perturbation (see Ref. [35], where a similar 



3 Moreover, we demonstrated in Ref. [7], that the deformed parafermionic currents, acting 
on the deformed vertex operators, reproduce also VEVs for perturbed W algebra primaries, 
and, in particular, explain the factorized form of resulting VEVs H56II of the fields Q^^- 
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phenomena was studied in the sine-Gordon model context). It can be directly 
checked, that for perturbed fields we have relations 

A m n k = , A\ +k Mfe = . (38) 

N -JJ- 

The first equation here follows from the parafermionic currents action procedure, 
introduced in Ref. [7J. Indeed, one can check, that taking the a — > limit in 
the expression 

<A^$i%J{/3},{/3'})(n-i,«)~ (39) 

Mme-^««n B C9i)?-*II Bt Cfi5)2-* B (°0 ' 

we get exact zero for arbitrary particle number n. For the second equation, 
after the application of the Wick theorem, we find that the condition 

(QI11+, /i fe |/3 1 ,..., ( s+1 ,/3i,...,^)=O, (40) 

N 

for the arbitrary 2s + 1 particle matrix element of the scaling field A\ +k fi^ is 
reduced to proving the following trigonometric identity 

e ff * n h e -- (fe+i)6 ' n n «**(^ - + w * 

{«y },{&.,} i=l J=l i=l j=l 

smh( ^ + — (a* - o J -)ll«nh( ^ " - &i) - 0, 

which follows from equation (|2"51) . 

Still, the null vector conditions, as well as other relations in the the space 
of form factors of scaling fields, have to be studied deeper, due to a possible 
appearance of the fields, satisfying the conditions ©-(HI)- The following example 
illustrates, that we have to be very careful with an analysis of the off critical 
fields structure. We consider the multi particle form factors of the perturbed 
fields 

tPk<P k — v 2k-N,N-2k • 

The expressions for these form factors, computed by a direct application of eq. 
(|28[) . up to a normalization, coincide, as functions of rapidities, with the corre- 
sponding results for the field Q^k -2k = ®2k-N N-2k- The equality follows from 
the identities for trigonometric functions. In that way, after the perturbation, 
we get, formally, that two fields have the same form factors, while they are dif- 
ferent at the criticality. This disagrees with the statement, that conformal and 
massive scaling fields should be in one-to-one correspondence [13]. The point 
is that, this formal coincideness of expressions for different fields happens in a 
very specific case, in which the form factor prescription has to be worked out 
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more carefully. Using the explicit values for the conformal dimensions (pH) . we 
see, that the following relations between the scaling dimensions of these fields 
takes place 

24* - 2d 2k =(2 -2 £>i). (41) 

Since Dj = A £l is the conformal dimension of the thermal operator e\, this 
equation is exactly the condition , that the fields ijjki> k are in the first order 
resonanc^l with fields ^k -2k- ^ there is an ambiguity (0} in its definition 

^ki>l ~ + const A $ { 2k,-2k ■ ( 42 ) 

Moreover, the spinless fields tpk'ipk formally have divergent expectation values, 
therefore, we need to work out a regularized prescription for its normalized multi 
particle matrix elements. 

To introduce a well-defined form factor of the field V'fcV'L we propose to use 
a freedom in the choosing the index N. Namely, let us provide an analytic 
continuation of the expression for the multi-point form factor of this field, by 
changing ./V — > N + e, where the regularization parameter e is a small number, 
which will be set to zero in the end. After the change, the vacuum expectation 
value becomes finite and proportional to 1/e. We can divide the e-deformcd 
expressions, obtained from eq. ([251) . to the correspondent finite value of VEV 

and normalize the null point form factors to 1. The same analytic continuation 

(2) 

and normalization are provided for the field & 2 k _2k- Then, we propose to define 



the regularized form factors of V'feV'fc as a result of the formal calculation 

(43) 



lim — 

e->o e 



V <^> w <*&U> (e) 



This is one of the possible prescriptions in fixing the ambiguity (j4|) . In the ultra 
violet regime, we should provide an analogous regularization for the field V'fcV'l-i 
which leads to vanishing expectation values of regularized fields. 

To support this construction, we applied the regularization procedure for 

(3) ~ (3) 

findingthe form factors of the descendant field E\ = W_{W_{s\ in the N — 4 
model|£| The field E\ in this theory satisfies the first order resonance condition 
with the identity operator. For form factors of this field, the regularization 
prescription provides a correct result, which can be calculated from another 
consideration. Namely, for general N, the form factors of the field E\ have a 
very simple property. It can be expressed in terms of form factors of the first 



4 The same is true for the Z?j charged fields V'feV'fe an( l ^2h2h' 

5 This theory can also be considered as a well-known sine-Gordon model in the reflection- 
less point, with the parameter /3 2 = i . The form factors of the exponential fields in this case 
can be extracted from the results of Refs. 1221 1331 
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energy field in a following way 

m ' wlM ° s s (■-»+«-*) x 

x( £l |{/3},{/3'}) (liri) . ' 7 (44) 

This equation comes from the fact, that in CFT [27] the energy field e\ has a 
null vector at the level 2. This means, that there is a linear relation between 
the descendants [TH] W_ 2 ei and L-iW_{e\. Now we have to take into account, 

that the modes W_ 2 of the W^z) current acts on the energy field E\ as the spin 
2 integral of motion. Due to its even spin, it is odd with respect to the charge 
conjugation transformation, and have the form 

n 

(W^exm, {/3'}> (n , n) = E - e2 "0 (^I^M/^W) ■ ( 45 ) 

This completes the explanation of the cq. (|44|) . 

We suppose to study another relations, following from the null vector con- 
ditions, in a separate publication. 



4.3 Explicit expressions for form factors 

Finally, let us derive explicit expressions for the form factors, which will be used 
in the correlation functions studies. With the definitions (|2"T1) - (|31I) . the equation 
(28) for the disorder fields has a conventional form 

E IT ajbje^> ( - ( - k+1 ^ aj ~ bj ^~ m( ' aj+bj ^ x 
{ a j, b j} i 

((Z ai (/? 1 )...Z b t J/3;j))(< ) _ m )- (46) 

As an example, we write down the explicit expressions for the first form factors 
of the disorder operator 

<A*fe]0> = , 

Gufclft^iht! = -C f (/3i - /3' 2 )(n k ){k}e^^-^ , 



(/i fe iA,/3 2 , m,i3' 2 ) iwn = c(/3i2)c(/?i 2 ) n - 



f { i /-, -i- — , , 

{!}{fc} {!}{fc} 



/ {2}{fc + l} {2}{fc-l} \ 



1G 



where Vj and Tj are j-th symmetric polynomials of the variables and e^, 
respectively, and we introduced the notation {a} — sin(^). In eqs. (l47l) . we 
have written explicitly the normalization of scaling fields to its VEV [33J . The 
exact values of VEVs for the relevant for us fields [7 will be written below. It 
is rather direct task to find out expressions for the higher particle form factors. 
However, in our numerical computations, we will not need it. 

Another simplest example is the case of parafermionic currents. Our candi- 
date for the parafermionic current form factors can be easily found by using the 
equation (|28f 

2 

Wi,#,#)itii = -^{2}C(/3 12 )nC t (/3i -#)x 

3 

] ' 3 

W/9i,A,^,^,^ 3 >ititm = i ^{2} 2 C(/3i 2 )nC(^-) 

i<j 

2 3 

x nflC f (ft - fyeG-iM+M-W+K+M 

i 3 

x(r 2 + i/ir 1 + ^ 2 (l + 2cos^)) , (48) 

etc. By we denoted the normalization factor, determining the one particle 
form factor. We will fix it explicitly in eq. (|7Tj) . 



5 Correlation functions 

In this section, we develop the conformal perturbation theory and compare long 
and short distance asymptotics of correlation functions of scaling fields. Our 
aim is to demonstrate, that the asymptotics are in agreement with each others 
at the intermediate distances, and, therefore, give an effective description of 
correlation functions at all scales. See Refs. [H HOI [3TJ [H] for other results 
in this direction. 

Before proceed further, let us make an important comment. The long dis- 
tance behaviors of correlation functions are described in terms of the mass of 
the lightest particle parameter. The exact relation between the mass and the 
coupling constant can be derived via TBA technique [2]. For the sine-Gordon 
model, this was done in the Al. Zamolodchikov's paper [3J. In our case the 
mass-coupling constant relation can be found from the results of Ref . [35] . To 
simplify expressions, here and below we use the notations: 

r(a) i r(2/A)r(i - i/n) , , 

7 (o) = , u = , k = M , , ' ■ (49) 

n ' r(l-a) N + 2 ' r(l/A0 v ' 

Applying these definitions, the explicit relation between the parameters M and 
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A is given as following 

(2ttA) 2 = K^- 2u ^(ti)j(3u) . (50) 

5.1 Conformal Perturbation Theory and exact VEVs 

In this subsection, we consider the basic notions of the conformal perturbation 
theory [I] and give explicit values of the VEVs [3 [8] of the fields, which we will 
need in the computations. 

In what follows, we develop the conformal perturbation theory pQ for the 
two point functions of scaling fields $ a and $b 

(* o (M)* 6 (0))=£c£i(;M) (0,(0)), (51) 



where C a l b (z, z) are the structure functions, and operators Oi form the basis in 
the space of scaling fields. Functions C^(z, z) can be expanded in the pertur- 
bation series 

C° l b {z,z) = | z | 2A o lz -A* a -A„ b ^-A* a -A^ x 



where the coefficients C^l° l are the structure constants from the conformal 

fields theory, while the first order corrections C^l° l can be expressed through 
the integrals of correlation functions in CFT 

GiT = J d "y (*a(O)«Ml)ei(0, v)Oi(«>)> ■ (52) 

The vacuum expectation values (O/) of the scaling fields Oi, appearing in the 
equation (f5"Tj). have a non-perturbative nature pQ. These fundamental quan- 
tities depend on the normalization prescription for the fields. We find VEVs, 
assuming that the scaling fields are satisfy the standard conformal normalization 
prescription 

(0(2,2)6^(0)) = |z|- 4Ao +••• , \z\^0. (53) 

The exact VEVs for physically important operators ^>m]m an d -E* = wf^ W 

in scaling Zm Ising model were found in Refs. [3|8], following the approach of 

Refs. [4, 35, 26] . The vacuum expectation values for the disorder fields are 

(/i fc ) = n 2d » x (54) 

f°° dt ( sinh hut sinh(iV — k)ut , +. 

exp / — ^ 2d k e M 

J t \ sinh t taah N ut 
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It is interesting, that for even values of k, the integral above can be calculated, 
giving 



For a general case, _ m = (ip^) 1 (ip) 1 fik, where m = k — 21, the vacuum 
expectation values of fields are elegantly expressed in terms of (/%) as 



(*ff-J _ / K (iV + 2)\^^( l + l)7( !±i 



</*} V * y ii (fc -z) 7 (^i) 

In particular, VEVs of thermal fields e k = (ip) k {ip) k ^2k correspond to the m — 
case in the equation above: 

^(wywrmgrn'* (57 , 



TV / (2fc)! V 7(«) 



7 2 (^r) 



S"7(^)7((2* + l)u) 



We also obtained the exact result for the expectation values of the normalized 
descendent fields E k = W^W^le k 

(E k ) 2 (iv + 2) 2 r 2 (i + ^i)r 2 (i-f) 



(e*) ~~ ' 27V r 2 (l-^±l)r2(l + A) 

r(i + f)r(i-^±g) 
r(2-f)r(2 + ^± 2 )' 



(58) 



5.2 Correlation functions (/ii(x)//i(0)) 



In our previous paper [7] , we studied long and short distance asymptotics of the 
correlation functions (ai(x)a\(0)) and (//i(x)//|(0)) of order and disorder fields. 
For completeness, we collect the correspondent ultra violet expansion data in 
the appendix A. We found, that IR and UV asymptotics of the correlators match 
at the intermediate distances. This confirm our expressions for VEVs, as well 
as form factor expressions. 

In this subsection, we want to discuss a more complicate case of the cor- 
relation function, including two disorder fields (/ii(x)/xi(0)). We found, that 
studying this case is rather instructive, since it gives an example, where the 
resonance fields ©-© appear at the short distance expansion for all integer 
N = 2,3,... 

The long distance expansion of this correlation function is provided in a 
standard way. In a two particle form factor approximation, we have the spectral 
decomposition 
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+{ t i(z,z)\9 1 ,h) 1 t il (M + (0)|e3+t7r,di+i7r)it,iJ +••• . (59) 

Studying the short distance behaviors is more involved. Indeed, we have the 
following leading terms in the conformal perturbation theory expansion (|51l) : 

(|*i(;M)/il(0)) = C^ 1 (r)(^)+C^ 1 (r)(r ] ) + 

+Ctt 1 {r){^)+--- , \z\=r. (60) 

(2) 

Where we introduced the notation r\ for the field r\ = $ 4 2 ■ The leading contri- 
bution to the correlator (|60l) at the short distances comes from the zero order 
term C^^ 2 (r)(/j 2 ) of the perturbation theory. It can be computed using CFT 
structure constants, found in Ref. [9], and the exact vacuum expectation value 
for the second disorder field /12 (|54|) . 

Cj»£{r)M = r-*^(rK)^(^y , „ = ^ . (61) 

At the first order perturbation theory we have the term, including the field 77. 
The contribution from this field to the perturbative expansion can be expressed 
in terms of the product of the CFT structure constants Cffi^ Cj^ei , multiplied 
by the simple two dimensional integral over y from the correlation function 

( W (0) ft (l) El (M)r,(»)) CfT = M 2 "|l - y\ 2u . 

Using the exact result (1551) for VEV of the field 77 = $4 -2> we obtain the analytic 
result for this term in the first order approximation 

ig [tk) (l + 2u )2 7 (3 U ) 7 2 (2 U ) [nu) ■ (b2) 

A more difficult task is to compute the first order correction to the structure 
function Cfff ( r ) anc | the contribution, coming from the term ipi/)K We note 
that the conformal dimensions of the fields ipffi and \i 2 satisfy to the relation 

N 

2A V ,-2A M2 =2-^=2(1-A £l ). (63) 

The first order resonance condition ([3]) appears now at the short distance ex- 
pansion. The equation (|63p basically means, that two terms in the perturbation 
theory will have the same powers of r. Whenever this happens, we expect, 
that divergencies coming from the contributions from the term C^f^r)^^) 



20 



should be canceled by divergencies, appearing from the first order correction to 
the structure function Cft)jt^{r){fJb2)- This phenomena, in general, should lead 
to logarithmic terms in a perturbative expansion. 

The appropriate contributions to the resonance terms can be computed 
within the Zjq Ising model settings by providing the analytic continuation in the 
parameter N, as it was discussed before. I.e., we can compute the contributions 
after the change N — > N + e, and then consider the limit e — > 0. To check 
the validity of this approach, we also perform the computations in a different 
method. We use the fact, that the parafermionic CFT with the central charge 
((T|) belongs to the series of W algebra symmetric unitary minimal CFT WM 



(p) 

N 1 



where the parameter p is chosen to be p = N + 1 [37]. The field theory @, 

from the view point of this model, is the simplest representative of the series of 

(p) 

integrable perturbations of WM£ minimal models by the primary field with 
the conformal dimension = 1 — which, in our case, coincides with the 
first energy operator. To study contributions of the resonance fields we can do 
the deformation of the WMjj' by changing the parameter p in such a way that 
the variable u=l/(p+l) would have the form 

1 



N + 2 

For small non-zero e the expressions for the resonance fields contributions are fi- 
nite and well-defined. The result of computations can be written in the following 

way 

x 7 (j) 7 (Nue) Z ("') ^ ggj ^ {^T^) 7 (ggg*) 

x (1 + e(N + 2) X0 + 0(e 2 )) , (64) 

To simplify the resulting expression, we introduced a shorthand notation for the 
term in the last line 

(N~2)u\ . (\ (N-2)u\ 2 



2V \2 2 ) \2 2 / JV 

'N + 2' 



+ (iV + 2)log 



N 



The symbol is reserved for the Euler's constant and if> stands for the loga- 
rithmic derivative of the gamma function ip(x) = T'(x)/T(x). (Unfortunately, 
this standard notation conflicts with the usual symbol, which we choose for the 
parafermionic current ip(z)). 

The contribution, which comes from the term proportional to (i/jffi), is found 
to be of a similar form 



2 

X 
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/ 1 y"" 1 -f(Nue) ( \ - Nu € \ (l-{N-2)u t 

/ 1 i 1 \ [N + 2 1 
x 7 tt 1 - T7 7 



Nu e N J \ N Nu ( 

2> 



where 



x(l + e(N + 2) X i + 0(e 2 )) 1 (65) 



. , ,T Nu\ , fl Nu\ , 

Xi = -=N»h; + -s--^(n--s- + • 



Now we look for the regularized expressions in the limit e — > 0. 

It can be derived from the equations (|64|) -(|65 |) . that both terms have a single 
pole at e = 0. However, as we expected from a general consideration, the sum of 
residues at this pole is zero, and the total contribution of these resonance terms 
to the correlation functions is well defined in the e — > limit. Computing the 
limit, we find, that the correctly defined sum of two resonance terms, is given 
as 

x ^8( 7B - 1) + Wlogrnei™- 1 + 2(N + 2) ^(1) + 1>(~) 

- N («im )+ «w&>))- (66) 

The numerical data for long and short distance asymptotic expansions are shown 
in the Fig. 5 for the case N=7. 

M" 
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Fig. 5. Correlation function of disorder fields at N=7. 

The dashed curve is for UV asymptotics, while the full line denotes the form 
factor decomposition, up to two particles. We observe, that there is a matching 
of the asymptotics at the intermediate distances. We found, that, in the region 
0.01 < Mr < 1, the long and short distance asymptotics agree with the relative 
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error around 1 percent. This, rather good, numerical preciseness confirms our 
hypothesis on the identification of form factors, as well as the short distance 
regularization prescription. 

For other values of the parameter TV, we provided similar numerical computa- 
tions. We found, that the relative error decreases with increasing of the number 
N. For large N the error becomes smaller. For example, for N=ll case, it 
is already less than 0.1 percent. The agreement between data can be further 
improved by taking into account higher particles form factors, however this is 
beyond the scope of the present paper. 

In the limit of large parameters N we find, that our short distance expansion 
function behaves as 

(/il(z '^ l(0)) =1 + 4 f-2fl + (fi 2 9 ~)M*A +■-., (67) 
(Mi) iV V 2 J 

where f2 = log ^ Mre ^ B+ ^ . We checked, using the Mellin transform, that our 
form factor expression leads to the same expansion at the small distances up to 
the order iV~ 3 . In the next orders new corrections can appear from the higher 
particle form factors. 

Another non-trivial test for the correctness of our expressions is the limit 
to the Ising model point N — 2. In the case, when N — > 2, our ultra violet 
asymptotics leads to the following answer 

. . , , If Mr, Mre< B \ , , 

(MiMmi(O)) = -r (1 --J- log— — + ---J , (68) 

which agree with the known short distance expansion for the Ising model disor- 
der parameter correlation function |15i 136) . 

5.3 Correlation functions of parafermionic currents 

In this subsection we would like to compare short and long distance asymptotics 
for the two point correlation functions of parafermionic currents ip and ip^ at the 
region Mr ~ 1. We think, that a consideration of this case might be interesting. 
First of all, because the correlator of free fermions is one of the most simplest 
in the Ising model. Unlike the case of order and disorder fields, it is given 
exactly by the Bcsscl function, i.e. by the solution of a linear equation, which 
is simplier, than the Painleve equation [15, 36 . From the other side, this is one 
of few examples of correlation functions of operators with fractional spins. 

We recall, that in the unperturbed CFT, the currents ip(z) an d (V* ) have 
conformal dimensions (Ai,0) and (Ajv-i,0), defined by eq. (fP2"|) . In the ultra 
violet region the correlation function of these fields is expected to have a form 
(seeeq. (p]>) 

2)^(0)) = z~ 2Al +■■■ , |z|^0, (69) 

which is natural definition of the conformal normalization in the given case of 
operators with non-trivial spins. 
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5.3.1 General N case 



The leading contribution to the infra red asymptotics of the correlation function 
of parafcrmionic currents comes from the one particle form factor approxima- 
tion. We will see, that it already leads to the nice agreement between the 
asymptotics at the intermediate distances. Using the integral representation for 
the modified Bcssel function, we can express the one particle contribution in the 
following analytic form 

(#,z>t (0 )) = 2C^(r) + --, z = r. (70) 

For simplicity, we choose here and below the space coordinate to be zero, which 
corresponds to the real z = r. 

In the equation (|70l) we took into account, that for the fields with spin 
the vacuum expectation values vanishes. Our proposal for the exact value of 
the multiple , leading to the conformal normalization (|69|) of parafermionic 
currents, is 

r(l + £) (N + 2\ 2 ^ ^ 2tt 
C * ~ T(T—Lj {-IT ) S 2 (2n+-). (71) 

The constant S2(2n + ^) is defined by eq. in the appendix B. This ex- 

pression can be obtained, following the ideas of Ref. [33]. It comes naturally 
from analysis of the divergency in the VEV of the operator ipffi , which can be 
effectively provided in the framework of the iy-symmetric CFT, perturbed by 
adjoint field, by use of the results of the paper |35j . In another way, it can 
be obtained by analyzing the deformed parafermionic currents normalization, 
see Ref. [8]. Further, we will see, that this normalization coincides with the 
known exact results for N = 2 and N — > oo cases. We also establish numerical 
checks, by matching the long and short distance asymptotics for correlators for 
arbitrary N. 

From the general arguments of the conformal perturbation theory pQ, we 
found, that the short distance expansion of the correlator of parafermionic cur- 
rents has the form 

i 1 / 4jy 

(ip(z, z)^(0)) = -357(1 + (nr) 2 A 1 + («x)*r+2 A 2 + 

4(W + 3) 4 JV n 

+ (nr) N + 2 A 3 + (nr) «+ 2+ A 4 + 

4(jV + l) | n 6(JV + 4) \ 

+ (nr) «+2 +Z A 5 + (nr) A 6 H . (72) 

Here and below in this subsection, we use in the short distance expansions the 
notation r = \z\. For the computation of the coefficients Aj, (j = 1, . . . , 6) in 
this equation, it is convenient to use the quantum equations of motion (|3"Tf 

^ z ,z)=\^ ( 2 2 l(z,z). (73) 
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We recall, that the fields $20 = WO e i an d ^-2 = (' i /'^) e i! have different left 
and right conformal dimensions (A w ,A 6l ) defined by eqs. (0, (|23l) . 

To define the coefficient , we integrate over z the first term in the decom- 
position of the following two point correlation function 

(di>(z,z)^H<:,o) = x^{^ (z,z)^(<:,o). (74) 

In that prescription the coefficient Ai is computed in a simple way from the 
CFT three point correlation function 



($(>)^(l) £l(o o)) 



CFT 



multiplied to the expectation value of the field e\. An effective method of 
computation of other coefficients in the equation ()72|) is to integrate twice the 
series expansion for the following correlation function 

mz,z)W(C,Q) = ^<<«k^)$-?o(C,C)) ■ (75) 

More explicitly, starting from the conformal perturbation theory expansion for 
the two point correlator 

*)*- ( &(C, C)> = (J^X^ x (76) 

(<7 J (r) + C ei (r)<ex) + C^(r)(e 2 ) + C^{r){Bk) + C*(r)(e 3 ) + •••), 

we obtain, by integration over z and (, the coefficient A 2 from the coefficient C 1 . 
In a similar manner, the coefficients A 3 , A 5 and A 6 are related, correspondingly, 
with the zeroth order structure functions C £2 (e2), C El (Ei) and C e3 (£ 3 ). The 
coefficient A4 is related with the first order correction term C^ £l (si). 

Using the free field realization for conformal primary fields [HI E] , one can 
find an integral representation for the leading corrections to the structure func- 
tions in the expansion (fT5"j) . The integrals in the first order perturbation theory 
can be taken by applying the technique from Refs. [37] . The exact vacuum 
expectation values for the energy fields (efc), (k = 1,2,3) and E\ are given in 
cq. (1571) and (|55|) respectively. The results for numerical coefficients Ak can be 
written as 



^ 7 (f) ' 

(N + 2) 2 
A 2 = _L_1_; 7 ( u w 3w ) 

2 4N 2 (N - 2) ,V ' n ' 



A 



3 



(Ar + 2) 4 j(u) iHwhiw) 
"72(A^ + 3)(A^ + 4)Af 3 7 (3u) 7 (A) 
(^V + 2) 3 7 2 (f) 



^ = -16iV3(3iV + 2) 7 (£) 7(U)7(3U) JW ' 
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A 5 = - 



(iV + 2) 5 (iV-4) 1 7 2 (^) 7 (4 M ) 7 2 ( M ) 

8(N + l)(3N + 4) (N + 4) 2 (iV - 2) 2 N 7 (A) 7 2 (2m) 
1 (A + 2) 5 7 2(^) 7 2(^) 7 2(3 )< 



2400 NHN + 4)(N + 5)(N + 6) 2 7 (£)7(£M&) 
7 2 (2u) 7 » 
7 (3u) 7 (4m) 7 (5u) 



(77) 



Note, that the coefficient A4 was computed analytically only for particular val- 
ues of N. This is because the multiple Jn in this coefficient represents the 
contribution, which is given by the rather complicate integral over the plane, 
from the CFT correlation function 

Jn = ~/ d 2 y ($^(l)<i>_ ( 2 2 | (0) ei ( 2/ ,y) ei (oo)) CFT . (78) 

The asymptotics of the integral Jn, considered as a function of N, will be given 
for large values of N further, in eq. (|84p. 

In the next subsection, we discuss the particular cases of N, for which the 
correlation function have specific features, including resonances. 



5.3.2 Ising model N=2 case 

We consider, first, the consistency of our formulae for N=2 with the known 
results from the Ising model. According to the fusion rules in Ising model case, 
there is no higher energy fields in the short distance expansion. So, we have to 
restrict ourself to the two first terms coming with the coefficients A\ and A2. 
A formal substitution N = 2 leads to a divergency which is, of course, expected 
from our previous consideration of the resonance fields. When we consider the 
limit N — > 2, we find, that the residues at the pole are canceled, and the 
resulting expansion for correlator has a logarithmic scaling. Let us assume, that 
z = r is chosen to be real. Then, the N — > 2 limiting value for correlator will 
have the form 

(tp(z, z)ip(0)} = Urn r [I + [nr) A\ + («r) N + 2 A 2 H J = 

1 / (Mr) 2 , /Mr ^ i\ \ 

This equation agrees with the form factor ultra violet expansion, where the 
exact correlator is given in terms of the modified Bessel function K\(r). In 
the numerical computations, we can formally put N in the general coefficients 
Ai and A2 to be close to 2. Then, the long and short distance expressions 
would match at the intermediate distances, as it is expected. Since the infra 
red expansion result is exact one, it is instructive to see for our short distance 
formulae and efficiency of the analytic continuation prescription. For example, 
the curves, depicted at Fig. 6. corresponds to N=2. 00001. 
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0.25 0.5 0.75 1 1.25 1.5 1.75 

Fig. 6. Parafermionic correlators at N=2. 00001. 

In the figure the dashed lines is given by the first three terms in the UV decom- 
position. The full line, here and below, represents the one particle form factor 
expansion. The dependence on N, in a vicinity of the Ising model case value of 
N = 2, is smooth. Matching between long and short distance asymptotics of the 
Ising model case serves as one of the confirmation, supporting the normalization 
constant Cl (fflj). 



5.3.3 Three states Potts model N=3 case 

The model for N = 3 coincides with the scaling 3-state Potts model. Up to 
the normalization factors, the exact form factors for parafermionic currents (as 
well as for other primaries) were computed for this model in Ref. [29] . The 
correlation functions of order and disorder fields in this model were also studied 
inRefs. [HIS]. 

According to the fusion rules for N=3 case, there is no energy fields £2 and 
£3 and we have to omit terms with the coefficients A 37 A 6 . We show at Fig. 7, 
that even first three terms in the short distance expansion lead to a rather good 
agreement between correlation functions asymptotics. 

jyr 2 &i <w + > 




Mr 



Fig. 7. Parafermionic correlators at N=3. 
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In this case we can advance in the exact computations further. In particular, 
the coefficient A4 can be calculated exactly, since the integral J3 can be taken 
analytically as 

" 3 - 16^3(1)1^-1) • (80) 

Then, the contribution from this first order perturbation theory term is given 
by the expression 



7 2 (|)7 4 (i) 



N=3 , z = r 3. 11- 2« 7(1)76(1) 



(r/s)*. (81) 



It is possible also to compute the second order perturbation correction to the 
identity operator in the expansion (|72|). as well as the contribution of the descen- 
dant field TT, those VEVs can be computed from the results of Refs. [38l 139] - 
With that corrections, and with the term with the coefficient A§, taken into ac- 
count, the long and short distance expressions agree up to the distances Mr=3. 

5.3.4 N=4 model 

The model for N = 4 describes a particular case of the Ahskin- Teller model. 
For the parafermionic fields, proceeding as before, we find, that the resonances 
appear at the short distance expansion, as well. 

Namely, providing an analytic continuation of general expressions with pa- 
rameter N, we can easily see, that the term with the coefficient A3 diverges at 
N — > 4. However, the divergency, coming from this term, is canceled by the 
singular part of the term with the coefficient A4, since the integral Jn has the 
following expansion, when N approaches 4 

J N = I^L_ + l- 8+ 0(N-4). (82) 

In this case we meet the resonance condition © between the operators 62 and 
ei. As in the Ising model, let us look for an analytic continuation in N of the 
general expression. Finding the regular, in the parameter N, expansion, we 
obtain 



lim ((ftr) 4 ™+2y4 3 + (nr) 2+4: «+ 2 A4) = 

3 j(h) , .14 fTK -yE__27\ 

= 2^7 7^y (rK)3log lY e 28 ) ■ 



(83) 



Again, the presence in the expansion fields, which are in the resonance, leads to 
logarithmic dependence of the short distance approximation. Let us note, that 
there is no divergency at the term A§ 

lim ( K r) 2+4Nu+4u A, = 1 3 '' 35 T 2 (-)( K r)^ . 

' 2 3 -2 10 -5n 3 A ' 
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Substituting that regularized expressions into the short distance expansion, we 
get a well defined expansion. The numerical data (where the contribution from 
the second order correction was taken into account), are shown at Fig. 8. 

M- 2ai <W> 



0.5 1 1.5 2 

Fig. 8. Parafermionic correlators at N=4. 

The relative error between asymptotic values in the region 0.0001 < Mr < 1 in 
our numerical computations is less than 1 percent. 



5.3.5 Large N case 

The nice feature of the Zjv models is that the correlation functions in the ultra 
violet and infra red regions (|72[) are in agreement for an infinite set of models 
including those with an arbitrary integer N > 2. We propose in the large N 
limit the following approximation for the multiple J/v in eq. (1771) : 



T 18 

Jn ~ 5N 



1 
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(84) 



Taking the A4 term into account, we have, for the asymptotic expansion for 
correlation functions at N — > 00 in the vicinity of r — 0, the following expansion 



M" 2Al <^(M)r (0)} 



1 



1 



= 3 (T (Mr) + -T 1 (Mr) 



(85) 



Here, the functions Tq(x) and T\(x) are given by the expressions, which agrees 
with the one particle form factors formulae for large N and small scales 



T (x) 
Ti{x) 

+ 16 



1 



x 

T 



= 2 log X 
2 6og $ 



16 

x 2 

— (1 

2 V 



log x + -f E 
+ \ogx) + 



-J lo g: r + 2kV 2 



51og- 



13 

T 



(86) 



With that conventions, both, short and large, distance expansions produce sim- 
ilar looking curves for N > 6. The example of the correlation functions for the 
case N = 11 is given in the Fig. 9. 
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0.4 0.6 0.8 



Fig. 9. Parafermionic correlators at N=ll. 

To show the difference between the asymptotics, we draw also the re-scaled 
correlators r 2 ~ <*~ [tpij)') in the Fig. 10. 



2 ( M-l ) 

r » <i/ri/r> 




M r 



Fig. 10. Re-scaled parafermionic correlators at N=ll. 

Let us note, that the relative error for ultra-violet and infra-red asymptotics at 
the region 0.000001 < Mr < 1 become smaller with increasing of the number 
N. For example, for N = 7 model, the error in this region is 0.8 percent, while 
for N — 20 it is already 0.08 percent. 

To deal with that small N cases we considered the A4 term more accurately. 
We provided the Pade approximation for the integral Jm (l78l) between the known 
and fixed points N = 4 and N = 00. The details of the computations are 
collected in the appendix C. Using this, we found, that the matching between 
the asymptotics of the correlation functions for small N is within 1 per cent at 
the distances 0.000001 < Mr < 1. 

6 Concluding remarks 

In this paper we studied correlation functions of disorder fields and parafermi- 
onic currents for scaling Zn Ising models and found that the long and short 
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distance asymptotics approach each others at the intermediate distances. From 
one side, a matching between ultra violet and infra red asymptotics gives an 
effective way of studying basic behaviors of the correlation function of the theory 
at all distances. From the other side, it confirms our construction of the form 
factors of the scaling fields. 

We discussed algebraic relations in the space of form factors and outline 
the set of problems related with the form factors of the descendant fields. 
In particular, we stress the role of quantum equations of motion, by showing, 
that they appear naturally in the form factor prescription. From the other side, 
we demonstrated, that the equations of motion can be a very powerful tool in 
studying the ultra violet asymptotics, within the conformal perturbation theory 
Their application allows to find coefficients in the short distance asymptotics, 
which are rather difficult to study by direct methods. Another useful method, 
which we applied in the analysis of form factors and conformal perturbation 
theory, is the W extended symmetry of the model. The W symmetric models 
of CFT and their integrable perturbations are, in general, rather complicate 
and we found that it is interesting, that correlation functions can be effectively 
studied for such models, at least, in the particular cases. We also stress the 
necessity of more deeper understanding the fields, which satisfy the resonance 
condition. We have shown, that these fields have unusual properties in both, 
form factor approach, as well as in the conformal perturbation theory, and have 
to be analyzed carefully. 

In this paper, we did not consider ultra violet expansions for correlation 
functions of "heavy" fields O a , which are in resonance with some field Ob (J3|)-(jl|). 
If we use the proposed form factor regularization prescription for these fields, 
then the calculation of the short distance behaviors for correlators become more 
subtle problem, and terms, including log 2 r, will appear in the short distance 
expansion. The well-known example for this property is the correlation function 
(ipip^ip'ipfo)) in the Ising model. 

Another interesting phenomena in studying of resonance fields is that, some- 
times, the renormalization of primary fields is finite and does not require intro- 
ducing the parameter e for the correct definition of such deformed primary fields 
off criticality. We suppose to describe these problems in other publication. 
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A Short distance expansion for order and disor- 
der fields 

For completeness, we collect here also the results [7] about the following corre- 
lation functions 

G+(x) = (a 1 (x)a+(0)) , G_(x) = (m(x)i4(0)) , (87) 

of order and disorder fields in scaling Zm Ising model. We considered the leading 
terms in the short distance expansion 

G±{x) = r- 4dl (ci + C% + C £ ± 2 (e 2 ) + Cf 1 (E,) + ■ 



Up to the first order perturbation theory, the structure functions entering this 
expression are 

t _ A 7 rr 2 ( 1 - 2 ")( 7 ( M ) 7 (3 M )) 1 / 2 7 (4 M ) 
C ±- l± 2(l-4 U )2 7 3( 2u ) ' (89) 

r 4 "( 7 (u) 7 (3u)) 1 / 2 A7rr 2 w 2 7 (4u) 7 4 (u) 
C± — =F- 



2 7 (2u) 4 7 4 (2w) 
i2u A^r 2 ( 1 - 2 «)u 2 7 3 (u)( 7 (5u)) 1 / 2 
12(1 + 2m) 2 7 2 (2u)( 7 (3u)) 1 / 2 ' 



rEl _ 2+iu u(l-4u)( 1 (u) 1 (3u)) 1 / 2 



2(1 + 2u)(l - 2w) 7 (2u) 



Note, that the Ising model correlation functions short distance expansions [HI 
[39] agree with TV — > 2 expressions of these correlators as well0 

B Two particle form factors 

Let the contour C goes from infinity above the real axe, then around zero, and 
then to infinity, below the real axe. Introduce the notation 



S 2 (x) = exp- / — \ J log(-t) . (90) 



1 f dt sinh(x — 2ir)t 
i— / o 

2 J c 2-Kit sinh 2 (7rt) 

Then the functions appearing in the equation (I33|) . read 



i sinh(^) 
+ f ) sinh 



CCS) = 2 : ^-^5 — X (91) 

2sinh(| + f)sinh(|-f) 



S 2 2 (27r+f) 
C t (/3) _ _J ^(2,+ f) 



cosh I S 2 («/3 + 3^ + ^)5 2 (-i/3 + tt + ^) 



7 The same phenomena was observed by V. Belavin in the case 1401 of the minimal models 
of the CFT M2,2n+l perturbed first energy field. 
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C Integral Jjy 

In this Appendix we discuss briefly the problem of the calculation of the integral 
J N = I jG(z,z)d 2 z, (Jill), where 

G(z,z) = ($| (0) $ 2 _ 2j0 (1) e 1 (^,l)e 1 (oo)> CJ?T . (92) 

Using a free field representation for parafermionic CFT (see for example [SHE]), 
we can express this correlation function in terms of four dimensional integrals. 
It is convenient to apply the integral transformation derived in Refs. |45j and 
represent G(z,z) in the form 

G(z,z) = \z(l - z)\- 8u (Mi gi (z,z)-N 2 92(z,z)). (93) 

The numerical coefficients ATi and Af 2 are normalization factors. Their exact 
values are given as (we use the notation u = jfp^) 

N2 ~ n^(2uh(3u) ' M -^ 2 2(1-2 U )' (94) 

The functions gi(z,~z) and g 2 (z,z) in the equation (I93[) are defined by the 
integrals: 



9l (z,z) = / d%d 2 t 2 Q(t u t 2 ,z), g 2 (z,z) = I d 2 hd 2 t 2 R(t u t 2 ,z), (95) 
where 



R(t 1 ,t 2 ,z)^\t 1 t 2 ^-ti){l-t 2 )(z-t 1 )(z-t 2 )\ 2u \t 1 -t 2 \ iu 4 , (96) 

Q(h,t 2 , z) = R(ti,t 2 , z) ^ ■ (97) 

(<i - z)(t 2 - z) 

It follows from the results of paper [46], that the functions g\(z^z) and 52(2;, z) 
satisfy third order differential equations for each of the variables z and ~z. 
Namely, 

L\{z)gi{z,z) = L 1 (z)gi(z,z) = , 
L 2 {z)gi{z,z) = L 2 (z)g 2 (z,z) = . 

Here the differential operators L\{z) and L 2 {z) can be written as polynomials 
in the parameter u. For L\{z) we have the definition 

L x {z) = z(l- z)d 2 z(l- z)d z + u(8z(l- z)(2z-l)d 2 + l{ 1] d z )+ 

W{l (2) d z + 8(1 - 2z)) + u 3 48(l - 2z) , (98) 

with 
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l[V = 2(z 2 + (l-z) 2 -llz(l-z)), 
4 2) = 12z 2 + 12(1 - zf - 522(1 - z) . 

The operator £2(2) is defined in a similar form as following 

L 2 (z) = d 2 z(l-z)d z z(l-z)d z + u(8z(l-z)(2z-l)d 2 + 4 1) &) + 

+m 2 (4 2) <9 2 + 8(1 - 2z)) + u 3 48(l - 2z) . (99) 

The coefficients ijp and are 

4 X) - -8(z 2 + (1 - z) 2 ) + 38z(l - z) , 

,(2) _ ,(2) 

'2 — n • 

Contrary to the case when correlation functions satisfy second order (hyperge- 
ometric) differential equation we do not know at present how to calculate in 
the analytic form the integrals associated with correlation functions satisfying 
third order equations. However, it is possible to calculate them numerically. 
The simplest way to do this is to develop a series expansion for the solutions 
convergent inside and outside a circle \z\ — 1. We can represent integral Jn in 
the form Jn — jffi + 4°°' > w h ere 

J^ = -f G(z,z)d 2 z, J^ = -[ G(z,z)d 2 z. (100) 
nJ\z\<i ' n J\z\>i 

Here we consider the integral jj^ . The calculation of integral jjy follows the 
same steps. 

We choose the basis of solutions to the eqs (|98I99[) yi,Wi in the form: 



Vi(z)=z vt to\l + J2$ ) (y)z*) , w i {z)=z v ^[l + Y^df{w)ziy (101) 

It follows from equations (|9"51 that v\(y) — v\(w) = 0, 1^2 (y) = ^{w) + 1 = 
1 + 2u, v${y) = v 3 (w) — 6it. All the coefficients dj (y) dj(w) can be derived 
from the corresponding differential equations. These coefficients determine the 
coefficients in the expansion of the "conformal blocks" 



Yi(z) = z- 4il (l - z)-^ Vl = z -^+Mv) fl + ^2df\Y)zA, 

\ .7=1 " ' 



(102) 



and similar coefficients d^ (W) in the expansion of z 4u (l — z) iu Wi. The value 
of the integral jjy can be expressed in terms of these coefficients as: 
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ji N = -It f C l {Y l {z)W l (z)~r l Y l {z)Y l {z))d 2 z= (103) 

V C ^ )(r) 4t.(,)-.(^)(^) ~ ndf{Y)df(Y y 

ti \tj j + l-Au + ^(y) 

Here the structure constants C% and numbers r, can be easily derived from eqs. 
(|93I95[) as the coefficients before the corresponding singularities at z — > in the 
correlation function G{z,~z). They are equal: 

1 7i(l - 2u) 7 2 (4 M ) 7 3 (u) _ 5 7 (5m) 7 (m) 

1 2' 2 (l + 2u)(l-4u) 7 4 (2u) 7 (3u)' 3 9 7 2 (3w) 1 J 

j 1 2u(l-2u) 2 

and n = 1, r 2 = - (1+2 ^ )(1 _i u) , r 3 = r 

For completeness we give here the expression for . To write the solutions 
outside the unit circle it is convenient to introduce the variable s = 1/z and to 
do the same substitution in the differential operators L\ and L 2 . The operator 
L\ is invariant under this substitution and operator Li transforms to another 
operator Li. The corresponding "conformal blocks" , Yi(s) = (1 — s) u yi(s) and 
Wi(s) = (1 — s)~ iM Wi(s) will have the following expansion 



Y^a) = 1 + £ 4° (*V ), W^s) = s Ui ^(l+Y,df ) (W)sA, (105) 

where the exponents f{(y) were given above and viiw) = v\(y), ^(w) = vi{y), 
v 3 (w) = v 3 (y) + 2. 

The integral can be expressed in terms of the coefficients dj (Y) and 
df{W) as 

-i c (t ^' tY) ^~' 7' 4 " (r)4 " (y) ) . .loo 

where numbers n were given above and q x = 3, q 2 = - (1 +2 M ) , 93 = m^Jjfl+g^) ■ 
We note that "dangerous" term with i = j = 1 in the sum (|106p does 
not contribute because (Y) = d^ (W) = and the first two terms for i = 3 

(91 — /0\ ' 

containing d_%(W) and d_{(W) also do not contribute because these coefficients 
vanish. The infinite sums (|103|106|) for small u converge rather fast. The terms 
in these sums decrease as j~ a where a = min{4 — 8m, 3 + 4it}. 

The solutions yi, Wi and Wi as well as the corresponding differential operators 
admit the expansions in term of parameter u = jA^ ■ As an illustration we give 
here the expansion of functions yi, w\ and W\ up to 0(u ). 



35 



yi = 1 + 4ulog(l - z) + 12u 2 log 2 (l - z) + 

u 3 (241og 3 (l - z) - 16 log(l - z)PL(2, z) - 8PL (2) (2, z)) + 0(u 4 ), 

wi = l + 8wlog(l-z) + u 2 (281og 2 (l-z)-8PL(2,z)) + 

u 3 (i^ log 3 (l - z) - 481og(l - z)PL{2, z) - 16PL(3, z)) + 0(w 4 ), 
32 



184, , 32 



wi = 1 + 4ulog(l - z) + — ?r log (1 - z) 



« 3 (^r lo g 3 (! - *) ~ -ri^si 1 - *)pl% z) + Pi (2) (2, z))) + o( u 4 ). 



The functions PL(2,z), Pi(3,z), PL^&z) are Pi(2,z) = - /* f log(l - i), 
PL(3, z) = /; f PL(2,t) and Pi( 2 )(2,z) = / Q z f log 2 (l - t). 

Using u expansion for the solutions yi , Wi and Wi we were able to find ex- 
plicitly five first terms in the expansion of integral Jn 

J N = + 2u+ + a lU 3 + a 2 u A ) + 0(u 6 ) , (107) 

5 3 

where the coefficients a\ and a 2 have the form 

ai = ^(31^ 2 - 2(65 + 57C(3))), 

Q 

a 2 = --(234 - 32tt 2 + tt 4 - 67C(3) - 92<C(5)) . 

We can now construct the function J^ ade \ which has the expansion (|107[) . 
coincides with I 3 given by eq. (IBTJ1) at N = 3 and has the asymptotic ([521 at 
N ^ 4. 



T (pade) _ 4u 34m(1 - 1.032u + 14.312u 2 - 129. 08u 3 ) 

N ~ 3(1 - 6u) + 15(1 - 0.679w + 11.953m 2 - 78.82u 3 ) 



The function J^ ade ^ can be compared with the values J^ wm ^ , which were cal- 
culated numerically by Enrico Onofri |47] 



N 


T (pade) 
N 


j(num) 


3 


— 1 .2109 


-1.2109 


5 


1.6046 


1.6046 


7 


0.67566 


0.67522 


10 


0.4031 


0.4033 


20 


0.184705 


0.184721 


50 


0.0723001 


0.0723002 



(109) 



100 0.03603784 0.03603788 
As it is seen from the table, an approximation of the integral Jjy for N > 3 by 
the function J^ ade ^ is done with a very good accuracy. The maximum deviation 
between data is less then 0.06 per cents. 
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At the end of this Appendix we note, that it is possible to analyze in a 
similar spirit another integral, appearing in CFT perturbative calculations 

In= ^J d2z 2 \ z{i - z )i~ 8u U2 g2{z '^ > ( n °) 

(see eqs. (|93I95I) ')'). This expression appears as an integral of correlation function 
of four thermal fields ey. It determines the first perturbative correction to the 
structure constant C\\ 

In = iy"(ei(0) ei (l) ei (z) ei (oo)) CFT d 2 z. (Ill) 

The integral has an series expansion 
24 

In = —u(l + 2u + 8u 2 + byu 3 + b 2 u 4 ) + 0(u 6 ) , (112) 
5 

where the coefficients by and b 2 are given as 

by =2(31tt 2 - 220 - 54C(3)), 

b 2 = 2(196 - 34tt 2 + 2tt 4 - 67C(3) - 95C(5)) . 

The expression for j^ ade ' ) can be derived by using the same conditions as for 
calculation of the J^ ade ^ , besides the condition that /( pade ) coincides with I3 = 



10.0897. However the value j( pode ) = 10.0914 is very close to 7 3 . We obtain the 
following answer 

Apade) _ 8u 8u 8u(l + 13.053tt - 11m 2 + 107.2m 3 ) 

N ~ _ 3(l-6u) + (l-4u) 15(l + 1.053u- 11.585u 2 ) ' ( ' 

For N > 3 this function approximates Jjv calculated numerically [47] with an 
accuracy which is better then 0.03 per cents 



(114) 



N 


j-(pade) 
1 N 


r(nwm) 
*N 


5 


-0.2462 


-0.2462 


7 


0.56130 


0.56124 


10 


0.4640 


0.4640 


20 


0.240352 


0.240348 


50 


0.096099 


0.096099 



100 0.0480157 0.0480157 
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